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ABSTRACT 

HeH+ is found to be the dominant positive ion over a wide range of tempera- 
tures and densities relevant to helium rich white dwarfs. The inclusion of HeH+ 
in ionization equilibrium computations increases the abundance of free electrons 
by a significant factor. For temperatures below 8000 K, He~ free-free absorp- 
tion is increased by up to a factor of 5, by the inclusion of HeH"*". Illustrative 
model atmospheres and spectral energy distributions are computed, which show 
that HeH+ has a strong effect upon the density and pressure structure of he- 
lium rich white dwarfs with Te//<8000 K. The inclusion of HeH"*" significantly 
reddens spectral energy distributions and broad band color indices for models 
with Te//<5500 K. This has serious implications for existing model atmospheres, 
synthetic spectra and cooling curves for helium rich white dwarfs. 

Subject headings: stars: white dwarfs, stars: atmospheres, equation of state. 

1. Introduction 

Bergeron & Leggett (2002) analyzed the recently discovered white dwarfs SDSS J133739 
-1-000142 and LHS 3250 (Harris et al. 1999, 2001), identifying both objects as extreme helium 
rich cool white dwarfs. However they encountered significant problems when attempting to 
fit the spectral energy distributions (SEDs). Bergeron & Leggett (2002) concluded that 
the discrepancy between their SEDs and the observed fluxes, is due to the physics used to 
calculate their model atmospheres. Here we investigate the molecular ion HeH+ as part 
of the missing physics of helium rich white dwarfs. We demonstrate that the opacity of a 
helium rich white dwarfs photosphere is significantly affected by HeH+. From the discussion 
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of Fontaine et al. (2001), it follows that increased opacity arising from HeH+, will lengthen 
the coohng time for hehum rich white dwarfs with Tf.ff< 8000 K. 

The only attempts to study HeH+ in helium rich white dwarfs known to us was made by 
Gaur et al. (1988, 1991). They showed that HeH+ exists in significant quantities in hehum 
rich white dwarfs and suggested a search for the infrared lines of HeH+. 

2. Equation of state 

The equation of state (EoS) is a vital component of any model atmosphere, it links 
the state parameters such as temperature, pressure, density, and internal energy. It also 
calculates the relative abundance of each species within the gas, which are essential to 
obtain accurate radiative opacities. The photospheres of cool extremely helium rich white 
dwarfs have densities which can reach upward of 1 g cm~^, under such conditions the use of 
a non-ideal EoS is required. 

We have adapted the non- ideal H/He EoS of Luo (1997). This EoS accounts for the 
non-ideal effects of electron degeneracy, Coulomb coupling and pressure ionization, but lacks 
an accurate treatment of pressure dissociation. The abundance of H2 is estimated using an 
equilibrium constant for the reaction: H2 ^ 2H, so that H2 pressure dissociates as hydrogen 
pressure ionizes. 

To account for the pressure ionization of H~ we have added a term to the hydrogen 
ionization equilibrium, given by eq. (22) and (23) in Luo (1997), so that 

Vh- = Lh-/Lh 

Lh = Lhi + Lhii + Lh- (1) 

where uh- is the ionization fraction of atomic and ionic hydrogen in the form of H~, Lhi 
and Lhii are the grand partition functions of atomic hydrogen and a proton (see Luo (1997) 
eq. 23). The grand partition function of H~ is given by 

Lh- = Wh- exp(2A - Eh- /kT) (2) 

where A is the electron degeneracy. Eh- is the sum of the ionization potential of hydrogen 
and H~ (14.352 eV), and Wh- is given by eq. (11)-(16) in Luo (1997) using a characteristic 
radius for H~ of 1.15 A (Lenzuni & Saumon 1992). 

Under certain conditions, the trace ionic molecules H^, H^, H^, HeH"'', and Hc^ are 
responsible for nearly all the free electrons in a H/He gas. We calculate equilibrium constants 
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for the formation of H2, H^, Hj, H3 , HeH+, and He^ from atomic H and He, H~, and 
free electrons with the Saha equation. Subject to conservation of charge, and of H and 
He nuclei, the equilibrium constants and ionization fractions are used to construct 3 non- 
linear simultaneous equations. These 3 equations are solved using a multi-variable Newton- 
Raphson technique. In this way the number densities for each species can be calculated for 
any given temperature, pressure, hydrogen fraction and value of A. The internal partition 
functions we use are detailed in Harris et al. (2004), for HeH+ we use the partition function of 
Engel et al. (2005). A converged value of A is found by iterating over a further conservation 
of charge equation: 

CeT3/2Fi/2(A-ecc/A;T) = AThh + iVHeii + 2iVHeiii - iVa- + iVn^ - iVn- 

+7Vh3+ + ^HeH+ + A^He+ (3) 

where Nj. is the number density of species x, ecc is the free electron Coulomb coupling 
energy (Luo 1997), -F1/2 is a Fermi-Dirac integral, T is temperature and a constant Ce = 
(2^/^/7r^)(A;me//l^)^/^. The left hand side of eq 3 is the number density of free electrons (see 
Luo (1994, 1997)) and the right hand side counts the charge on all ions. 

Figure 1 shows the number fraction of the species within our EoS, as a function of H to 
He number ratio, density, and temperature. At 5000 K and density of 0.2 g cm~^ Hjj" is the 
dominant positive ion for the hydrogen rich case. HeH+ is the dominant positive ion for the 
helium rich range —10 < \og^Q{NH/NH^ < —2.5 and HeJ becomes the dominant positive 
ion for Xog-^Q^NH/NHe) < —10. Figure 1 indicates that HeH"*" continues to be the dominant 
positive ion over a range of densities and temperatures. Lenzuni et al. (1991) present an 
EoS and mean opacities for a H/He gas of 72% hydrogen by mass, they correctly state that 
the opacity coefficent of HeH+ is wholely irrelevant. However, as illustrated below, for a 
hehum rich mix HeH+ strongly affects the opacity and cannot be neglected. 



3. Opacity function 

The opacity of a gas under the extreme pressures found in the photospheres of helium 
rich white dwarfs remains in question (Iglesias et al. 2002; Bergeron & Leggett 2002). The 
opacity of a cool helium rich atmosphere is dominated by H2-He collision induced and He~ 
free-free absorption, and He I Rayleigh scattering (Malo et al. 1999; Iglesias et al. 2002; 
Rohrmann et al. 2002). As such the opacity is strongly dependent upon the abundance of 
free electrons and H2. The sources of opacity data that we use is discussed in Harris et al. 
(2004). 

The monochromatic absorption coefficient at p = 0.5 g cm~^, logio(-/VH/A^He)=~5, over a 
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range of temperatures, computed both including and neglecting HeH+ from our EoS, is shown 
in figure 2. It is evident that if HeH"*" is neglected the gas opacity can be underestimated by 
as much as a factor of 5, over a significant range of temperatures. The dominant opacity, 
across the frequency range shown in figure 2 and for temperatures upward of 5000 K, is He~ 
free-free absorption. At lower temperatures collision induced absorption, in the infrared, and 
He I Rayleigh scattering, in the visible/ultra- violet, become important and eventually take 
over from He~ free-free. 

To determine if HeH+ rotation- vibration lines would be observable in a helium rich white 
dwarf wc have employed the recent publicly available HeH"*" linclist of Engel et al. (2005). 
We find that the absorption lines of HeH+ are too weak to overcome the continuous opacity, 
under the temperatures and densities found in helium rich white dwarfs. Therefore HeH+ 
lines will not be visible in the spectra of helium rich white dwarfs. For a discussion of HeH"*" 
line opacity and some of the temperatures densities in which it is important see Engel et al. 
(2005). 

4. Model atmospheres &; spectral energy distributions. 

We use the plane parallel model atmosphere code MARCS (Gustafsson et al. 1975), mod- 
ified for the new non-ideal EoS subroutines, discussed in section 2, and the new continuous 
opacity subroutines, discussed in section 3. The new EoS and opacity function subroutines 
are fast enough to be run in real time. 

As discussed in Saumon et al. (1994) and Bergeron et al. (1995), in the optically thin 
regions, the unusual opacity function of a metal free H/He gas results in multiple roots in the 
equation of radiative equilibrium. The high temperature solution to radiative equilibrium 
in the optically thin regions is preferentially found in our models. Such a solution is not 
physically realistic, rendering our models of Te//<5000 K below logr^ = —2 unreliable. 
However, as this only occurs at very small optical depths the emergent flux is unaffected. 

We also experienced a problem with convergence of the convective flux at temperatures 
of 5000 K and below. The pressure-temperature gradient (V) is very close to the adiabatic 
gradient (Va^), so that (V — Vad)/V ~ 10^^ in the convective zone. In the cool highly 
non-ideal regions, numerical noise in the value of Va^ calculated within our EoS is of this 
order, resulting in convergence problems with the convective flux. We have therefore not 
been able to obtain converged models below Te//=4500 K. 

We have computed a set of model atmospheres for \ogg = 8, log]^Q(A^H/AHc)= 10"^, and 
between effective temperatures of 4500 and 8000 K, including and neglecting HeH"*". Figure 
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3 shows optical depth verses temperature and density for model atmospheres of 4500, 5000, 
6000, 7000 and 8000 K. Although the temperatures remain relatively unperturbed by the 
inclusion of HcH+, there is a very strong affect upon the density and pressure. If HeH"*" is 
neglected then the density and pressure can be overestimated by up to a factor 5, similarly 
the electron pressure significantly underestimated. For Tg// of > 8000 K there are significant 
electrons released from H II and He II, which reduces the importance of HeH"*". 

Figure 4 shows the SEDs of our 4500, 5000 and 6000 K models, with and without HeH+. 
The 4500 and 5000 K SEDs show a significant changes if HeII+ is included in the ionization 
equilibrium, but the effect is only small for the 6000 K model. The reason for this is that 
above ~5000 K He^ free-free is the only significant source of opacity, so although the total 
opacity is increased the shape of the absorption function and hence SED is unchanged. For 
temperatures below 5500 K, He Rayleigh scattering and He-H2 coUision induced absorption 
contribute to opacity. As these opacity sources are unaffected by the increased abundance 
of electrons from HeH"*", the increase in He~ free- free opacity changes the shape of the total 
opacity function and SED. These differences are reffected in the broadband color indexes 
given in tables 1 and 2. These colors were computed by using the bandpasses given by 
Bessell & Brett (1988); Bessell (1990) and calibrating using a spectrum of Vega. There are 
significant differences, at Te//=5500 K and below, between colors computed whilst including 
and neglecting HeH"*". The large increase in the V-K magnitude, and most of the other 
color indices indicates that the models calculated with HeH+ are significantly redder than 
the models calculated without HeH+, this is also apparent in the SEDs. In general all our 
colors are redder than the colors of Bergeron & Leggett (2002). 



5. Conclusion 

A non-ideal H/He equation of state (EoS) which includes the molecular ion HeH+ within 
the ionization equilibrium, has been presented. It has been demonstrated, that under helium 
rich conditions and over a range of temperatures and densities relevant to helium rich white 
dwarfs, HeH+ is the dominant positive ion. Using the EoS, we have computed a set of 
continuous opacities which illustrate that HeH"*" can indirectly increase the opacity of a 
helium rich gas by up to a factor of 5. Using the recent HeH"^ linelist of Engel et al. (2005), 
we have found that HeH"*" line opacity does not significantly contribute to the opacity at the 
densities found in helium rich white dwarfs. 

From a physical point of view one of the most interesting reasons for studying helium 
rich white dwarfs is that the densities of their photospheres access regions in which the gas 
is strongly non-ideal. Saumon & Chabrier (1991); Saumon et al. (1995) have studied the 
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pressure dissociation of H2 in a pure hydrogen environment. However, one of the shortcom- 
ings of our, and all other equations of state known to us is that there has been no study 
of the pressure dissociation of the important molecular ions, B.^. HeH+, and He^. Before 
we can fully understand helium rich white dwarfs, our understanding of the physics of cool 
dense H/He plasmas must be improved. 

Our EoS and opacity function have been incorporated into a version of MARCS (Gustafs- 
son et al. 1975). Using this code we have computed model atmospheres, spectral energy 
distributions and broad band color indices for an illustrative range of helium rich white 
dwarfs. We find that in all models below 8000 K the pressure and density of the model 
atmospheres is reduced by up to a factor of five by the inclusion of HcH^. Furthermore, 
HeH+ significantly reddens the SEDs and color indices, for models below Teff = 5500 K. The 
importance of HeH+ should prompt a review of all current model atmospheres, synthetic 
spectra and cooling curves for cool helium rich white dwarfs. 

We thank Prof. Bengt Gustafsson for providing us with a version of MARCS, Prof. 
Hugh Jones for providing a spectrum of Vega, and the UK Particle Physics and Astronomy 
Research Council (PPARC) for support. 
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Table 1: Color indices for models calculated whilst neglecting HeH+. logg' = 8 and 

logio(iVH/7VHe) = l(r^ 



2<// 


B V 


V R 


V K 


R I 


I J 


.J H 


H K 


4500 
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0.58 


0.48 
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-0.17 


5000 


0.72 


0.44 


0.81 
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0.35 


-0.18 


-0.22 


5500 


0.60 


0.38 


1.14 


0.36 


0.32 
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-0.04 


6000 


0.50 


0.32 


0.98 


0.30 


0.24 


0.12 


-0.01 


6500 


0.42 


0.27 
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0.25 


0.18 


0.09 
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7000 


0.36 


0.23 


0.59 


0.21 


0.12 


0.07 


-0.05 


7500 


0.30 


0.20 


0.43 


0.17 


0.07 


0.05 


-0.06 


8000 


0.25 


0.17 


0.30 


0.14 


0.03 


0.03 


-0.08 



Table 2: Color indices, for models calculated with HcH"'". log 9 = 8 and log]^Q(A^H/^He) = 10 ^. 
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